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First Detection of Ar-K Line Emission 
from the Cygnus Loop 
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Abstract 

We observed the Cygnus Loop with XMM-Newton (9 pointings) and Suzaku (32 pointings) between 
2002 and 2008. The total effective exposure time is 670.2 ks. By using all the available data, we intended 
to improve a signal-to-noise ratio of the spectrum. Accordingly, the accumulated spectra obtained by the 
XIS and the EPIC show some line features around 3 keV which are attributed to S He/3 and Ar Hea lines, 
respectively. Since the Cygnus Loop is an evolved (~10,000 yr) supernova remnant whose temperature is 
relatively low (<1 keV) compared with other young remnants, its spectrum is generally faint above 3.0 
keV, no emission lines such as Ar-K line have ever been detected. The detection of Ar-K line is the first 
time and we found that its abundance is significantly higher than that of the solar value; 9.0^3'g and 
8. 4+2 7 (i n un its of solar) estimated from the XIS and the EPIC spectra, respectively. We conclude that 
the Ar-K line is originated from the eject a of the Cygnus Loop. Follow-up X-ray observations to tightly 
constrain the abundances of Ar-rich ejecta will be useful to accurately estimate the progenitor's mass. 

Key words: ISM: abundances — ISM: individual (Cygnus Loop) — supernova remnants — X-rays: 
ISM 



1. Introduction 

The X-ray spectra from the supernova remnants (SNRs) 
show various emission lines of heavy elements blown off 
by the supernova (SN) explosions. Since these elements 
are generated by the nucleosynthesis process inside the 
progenitor stars, the metal abundance pattern provides us 
with a clue to obtaining the information about the type 
of the SN explosion and the condition of the progenitor 
star. 

The Cygnus Loop is one of the largest (~ 1°.4 in radius; 
Levenson et al. 1997) and the brightest SNRs in the X-ray 
sky. It is close to us (540tgo° P c 5 Blair et al - 2005 ) and 
the age is estimated to be ^10,000 yr. Although the SN 
type of the Cygnus Loop is still unclear, McCray & Snow 
(1979) proposed that the SN explosion had occurred in 
a preexisting cavity and some other studies also support 
this result (e.g., Hester & Cox 1986; Hester et al. 1994; 
Levenson et al. 1997). Hester et al. (1994) observed the 
Balmer-dominated northeast rim of the Loop and showed 
the blast wave was decelerated from ~400 km s _1 to less 
than 200 km s _1 in the last 1,000 yr, which suggests that 
the blast wave is now propagating into the cavity wall. 
Since the Cygnus Loop is an evolved shell-like SNR, the 
spectra from the rim regions should mainly consist of the 
shock-heated cavity material. Then, by using optical and 
X-ray data obtained from the whole rim region of the 



Loop, Levenson et al. (1998) calculated the size of the 
cavity wall and estimated the progenitor star to be B0, 
a ~15M star. This result supports the origin of the 
Cygnus Loop to be a core-collapse explosion rather than 
the Type la explosion. 

Meanwhile, based on the analysis of the ejecta emis- 
sion in X-rays, some previous studies also suggested that 
the Cygnus Loop is originated from a massive star of 12- 
15M (Tsunemi et al. 2007; Kimura et al. 2009; Uchida 
et al. 2009a). Tsunemi et al. (2007) observed this SNR 
along the diameter from the northeast to the southwest 
with the XMM-Newton satellite and showed the X-ray 
spectra consist of two components with different tem- 
peratures. They interpreted this result as follows: the 
low-temperature component originating from the cavity 
wall surrounds the high-temperature ejecta component. 
Comparing the metal abundances (Ne, Mg, Si, S, and Fe) 
of the ejecta component with the theoretical SN mod- 
els, they concluded that the Cygnus Loop is originated 
from the 15M star. Kimura et al. (2009) and Uchida 
et al. (2009a) observed other regions of the Loop and also 
inferred the progenitor mass to be 12-15M . These re- 
sults are consistent with the estimation by Levenson et al. 
(1998). Therefore the origin of the Cygnus Loop is most 
likely to be a core-collapse explosion of a 12-15M star. 
However, Tsunemi et al. (2007), Kimura et al. (2009), and 
Uchida et al. (2009a) show that the number density of Fe 
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(relative to that of O) is ~10 times higher than that of the 
core-collapse model. The other problem is that the central 
compact object such as a neutron star has not yet been de- 
tected. If the origin of the Cygnus Loop is a core-collapse 
explosion, these points still remain open questions. 

In this paper, we present the result of the spectral anal- 
ysis of the Cygnus Loop by using all the available data 
obtained by the Suzaku and the XMM-Newton satellites. 
The X-ray spectra show many emission lines, including the 
first detection of an Ar-K emission line from the Cygnus 
Loop. 

2. Observations 

The Cygnus Loop has been observed many times using 
X-ray CCDs. We summarized the 41 observations in table 
1 (for Suzaku) and table 2 (for XMM-Newton). All the 
data were taken between 2002 and 2008. Their fields of 
view (FOV) are shown in figure 1. Top and bottom panels 
show the FOV of the Suzaku XIS (Koyama et al. 2007) and 
the XMM-Newton MOS (Turner et al. 2001), respectively. 

All of the Suzaku data were analyzed with version 6.6.2 
of the HEAsoft tools. For the reduction of the Suzaku 
data, we used the version 9 of the Suzaku Software. The 
calibration database (CALDB) used was the one updated 
in May 2009. We used the revision 2.2 of the cleaned event 
data and combined the 3x3 and 5x5 event files. The data 
obtained after November 2007 (Obs. ID: 501017010) were 
all taken by using the spaced row charge injection (SCI) 
method (Prigozhin et al. 2008) that reduces the effect of 
radiation damage of the XIS and recovers the energy reso- 
lution. In order to exclude the background flare events, we 
obtained the good time intervals (GTIs) by including only 
times at which the count rates are within ±3<j of the mean 
count rates. We examined both a local background and a 
blank-sky data for the background subtraction. The local 
background spectrum is extracted from the entire source- 
free regions on the periphery of the Loop while for the 
blank-sky spectrum, we accumulated all of the Lockman 
Hole data obtained since the launch of Suzaku. We con- 
firmed that there is no significant difference between two 
results. However, the total count in the Lockman Hole 
is an order of magnitude larger than that in the local 
source-free region. We therefore prefer the Lockman Hole 
background to the local background. We note that the 
position of the Cygnus Loop (I = 74°, b = —8°. 6) appears 
to be beyond the Galactic Ridge X-ray Emission (GRXE; 
Sugizaki et al. 2001; Revnivtsev et al. 2006). The solar 
wind charge exchange (SWCX) is also considered to be one 
of the causes of the soft X-ray background below 1 keV 
(Fujimoto et al. 2007). However, in terms of the Cygnus 
Loop, the SWCX is also negligible because of the promi- 
nent surface brightness in the low-energy band. Thus, for 
the background subtraction of the Suzaku data, we ap- 
plied the Lockman Hole data. These observation dates 
were close to those of the Cygnus Loop observations and 
we confirmed that they have no background flares. 

All of the XMM-Newton data were processed with ver- 
sion 7.1.0 of the XMM Science Analysis System (SAS). 
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Fig. 1. ROSAT HRI images of the Cygnus Loop (0.1-2.0 
keV). The FOV of Suzaku XIS and XMM-Newton MOS are 
shown with solid rectangles (top) and circles (bottom), re- 
spectively. In top panel, red, blue, and black correspond to 
Region-A, Region-B, and Region-C, respectively. 

The current calibration files (CCFs) used were the one 
updated on June 2008. We used the data obtained with 
the EPIC MOS and pn cameras. These data were taken by 
using the medium filters and the prime full- window mode. 
We selected X-ray events corresponding to patterns 0-12 
and flag = for MOS 1 and 2, patterns 0-4 and flag = for 
pn, respectively. In order to exclude the background flare 
events of the XMM-Newton data, we determined the GTIs 
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Fig. 2. Spectra from the entire FOV of the Cygnus Loop. 
Top, middle, and bottom panels show spectra obtained from 
XIS, MOS, and pn, respectively. Background spectra are also 
shown in red. 

in the same way as the Suzaku data. After filtering the 
data, they were vignetting-corrected using the SAS task 
evigweight. For the background subtraction, we employed 
a blank-sky observations prepared by Read & Ponman 
(2003) for the similar reason with the Suzaku' s case. 

3. Spectral Analysis 

3.1. Line Features at 3.0 keV 

As shown in figure 1, our FOV covers more than half 
the region of the Cygnus Loop. We improved a signal-to- 
noise ratio and created a global spectrum of the Cygnus 
Loop by using all the available data listed in tables 1 and 
2. We should note that the exposure time is different 



from region to region and some parts have not yet been 
observed. Therefore, the spectrum compiled in this way 
is not the precise one of the entire Cygnus Loop. 

First, we created a spectrum from each observing re- 
gion and combined them with FTOOLS mathpha. In this 
analysis, we treated the XIS, MOS, and pn data indepen- 
dently. In order to generate response matrix files (RMFs) 
and ancillary response files (ARFs), we employed xisrm- 
fgen and xissimarfgen for the Suzaku data (Ishisaki et 
al. 2007), rmfgen and arfgen for the XMM-Newton data. 
We added the RMFs and the ARFs with FTOOLS addrmf 
and addarf, respectively. Thus, we obtained the combined 
XIS, MOS, and pn spectra. Figure 2 shows the spectra 
from the entire FOV of the Cygnus Loop. Left, middle, 
and right panels show spectra obtained from the XIS, 
MOS, and pn, respectively. We also show background 
spectra in red lines. In these spectra, we can see many 
emission lines: C Lya, N Hea, O Hea, O Lya, the Fe-L 
complex, Ne Hea, Ne Lya, Mg Hea, Si Hea, Si He/3, and 
S Hea. Since these lines are all located below ~3.0 keV, 
they are clearly detected above the background level. It is 
difficult to detect line emissions above ~3.0 keV because 
the X-ray emission is buried in the background level; the 
Cygnus Loop is a middle-aged SNR, and the electron tem- 
perature is much lower (0.1-0.9 keV; Uchida et al. 2009a) 
than those of the younger SNRs. However, in figure 2, 
some line structures are seen around 3.0 keV above the 
background level due to the improvement of the statis- 
tics. These structures are likely S He/3 and Ar Hea lines 
which have not been detected by previous observations. 
They are seen in all three spectra, and it is particularly 
prominent in the XIS spectrum due to the low background 
level of the Suzaku satellite. 

Figure 3 upper two panels show 1.5-5.0 keV spectra ob- 
tained from the XIS data with different models. First 
we fitted the spectrum with bremsstrahlung plus five 
Gaussian components as shown in figure 3 left. In this 
model, the Gaussian components correspond to the emis- 
sion lines of Mg He/3, Si Hea, Si He/3, Si He7, and S Hea. 
The parameters of the electron temperature kT e and the 
emission measure (EM = J n e nndl, where n e and nn are 
the number densities of hydrogen and electrons and dl is 
the plasma depth) in the bremsstrahlung model were left 
free. In the Gaussian models, we varied the line center en- 
ergy and the normalization while we fixed the line width 
at zero. We fixed the absorption column density Ah to 
be 2.3xl0 20 cm -2 that is an averaged value based on the 
previous analysis of the Suzaku and the XMM-Newton 
(Uchida et al. 2009b). The best-fit parameters are shown 
in table 3. In figure 3 left, the spectrum shows an ex- 
cess from the model at about 3.0 keV. Thus we added 
two more Gaussian components to the model (figure 3 
right). Accordingly the reduced \ 2 value is improved and 
the F-test shows that the former model is rejected for 
a significance level of 99%. Furthermore, the line cen- 
ter energies of the extra components are derived to be 
2936±g? eV and 3128±39 eV, respectively. We also fit- 
ted the MOS and the pn spectra with the same model 
(figure 3 bottom). The result clearly shows an existence 
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of two Gaussian components at 2920^55 eV and 3132+29 
eV, respectively. Since these values are statistically con- 
sistent with the expected values for S He/3 (2884 eV) and 
Ar Hea (3104 eV) lines according to the Astrophysical 
Plasma Emission Code (APEC; Smith et al. 2001), we 
conclude that the excess structure around 3.0 keV is the 
S He/3 and Ar Hea lines. 

3.2. Abundance Measurement of Ar 

The X-ray spectrum of the Cygnus Loop is known to be 
expressed by two components; a high-temperature ejecta 
component and a low-temperature interstellar medium 
(ISM) component. Tsunemi et al. (2007) observed the 
Loop from northeast to southwest and reported that the 
high-/cT e component varies between 0.4 keV and 0.9 keV. 
Uchida et al. (2009b) observed the rim of the Loop and 
showed that the temperature of the low- kT e component 
ranges from 0.12 keV to 0.35 keV. They found that the 
temperature and the ionization timescale of each compo- 
nent are not uniform in the Loop. In order to measure 
the Ar abundance, we basically follow the two-component 
model with the NEIver 2.0 spectral code that includes K- 
shell emission line of Ar, but we introduced an augmented 
NEIver 2.0 (ver 2.0+) model (c.f., Badenes et al. 2006; 
Park et al. 2010) which includes various inner-shell lines. 
In addition, we prefer a VPSHOCK model (based on the 
NEIver 2.0+) to a simple VNEI because the VPSHOCK 
represents a more general case of the ionization timescale 
parameter than VNEI. For the following analysis, we also 
set the offset of the gain to be free by using gain fit in 
XSPEC for both detectors. We note that the best-fit off- 
set values are ~ — 5 eV which is within the energy uncer- 
tainty of the XIS (±0.2% measured with the 55 Fe calibra- 
tion sources 1 ) and EPIC (5 eV and 10 eV for MOS and 
pn, respectively 2 ). 

Consequently, we applied a two-component VPSHOCK 
model (NEIver 2.0+) : In the \ow-kT e component, we fixed 
the values of kT e and upper limit of ionization timescale 
T U pper (a product of the electron density and the elapsed 
time after the shock heating) to the averaged values based 
on the previous limb observations (Uchida et al. 2009b). 
We also fixed the metal abundances (in units of solar) to 
the values based on their results. The values are summa- 
rized in table 4. We fixed other elements (He, Ar, and 
Ca) to the solar values. We varied the abundances of O, 
Ne, Mg, Si, S, Ar, and Fe in the high-/cT e component, and 
set those of C and N equal to O, Ni equal to Fe. The 
other elements (He, and Ca) were fixed to their solar val- 
ues. Other parameters such as kT e and T U pp er were all 
free. We also set the absorption column density, TVh, to 
be free. 

As a preliminary step, we applied the two-component 
model for every full-band (0.3-5.0 keV) spectrum ob- 
tained from each FOV. Although this model fitted the 
spectra well, we could not constrain the abundance of 
Ar. Therefore, we next applied the same model to the 
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Fig. 3. The 1.5-5.0 keV spectrum of the Cygnus Loop with 
two model fits overlaid. We combined all XIS data shown in 
figure 1 top. They are fitted with bremsstrahlung plus five 
Gaussian components (left) and bremsstrahlung plus seven 
Gaussian components (right). The bottom figure shows the 
result of EPIC data. Red and black represent MOS and pn 
spectra, respectively. The residuals are shown in the lower 
panels. The best-fit parameters are shown in table 3. 



The Suzaku Technical Description: 

http: / / www. astro, isas. ac.jp/ suzaku / doc / suzaku_td / suzaku_td. html 

EPIC Calibration Status: http:/ /xmm2.esac.esa.int/external/xmm_sw_cal/calib/index.shtml 
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Fig. 4. Le/t: Spectra obtained from all XIS data. Right: 
Spectra obtained from all EPIC MOS data (red) and EPIC 
pn data (black). Each spectrum is fitted with the two-compo- 
nent VPSHOCK model. The residuals are shown in the lower 
panels. The best-fit parameters are shown in table 5. 

spatially-integrated spectra in order to measure the Ar 
abundance. We then found that the fits are far from ac- 
ceptable from a statistical point of view. The most notice- 
able discrepancies between the data and the best-fit model 
were found around 0.8 keV where the X-ray spectrum 
in figure 2 shows the highest count rate (which is about 
1000 times higher than that around 3 keV). As mentioned 
above, the low energy band (<1 keV) is dominated by 
the low temperature component, whereas the high energy 
band is dominated by the high temperature component. 
It is almost certain that the Ar lines come from the high 
temperature component and that the low energy band are 
not important to measure the Ar abundance. Thus, we 
here focus on the X-ray energy range above 1.3 keV. The 
discrepancies found at the low energy band is very inter- 
esting. However, detailed investigations are beyond the 
scope of this paper, and are left as our future work. As 
for the EPIC data, an excess exists around 1.5 keV which 
corresponds to the neutral Al-K line emission due to an 
instrumental fluorescence. Thus we used the data above 



Fig. 5. Best-fit parameters of Mg (triangle), Si (X), S 
(square), and Ar (circle) with the two-component VPSHOCK 
model are plotted as a function of O abundance. The values 
of x 2 are overlaid. 

1.6 keV for EPIC data. In these fits, we fixed the abun- 
dances for O, Ne, Mg (only for EPIC spectrum), and Fe 
to the averaged values of the results calculated from every 
FOV. 

Figure 4 left shows the fitting result of the XIS spec- 
trum. The best-fit parameters are shown in table 5. 
The spectrum is well fitted with the two-component 
VPSHOCK model (x 2 /dof=1096/1004). By comparing 
with several previous observations (e.g., Tsunemi et al. 
2007), it is reasonable to consider that the temperature 
of the high-/cT e component, 0.69^o'o2 keV, is an averaged 
one of the Cygnus Loop. The values of Mg, Si and S 
abundances are also consistent with the previous studies. 
From the results, it is clear that the Ar abundance of the 
Cygnus Loop is significantly higher than the solar value 
(9.0l3 g). We note that the calculated Ar abundance may 
significantly depend on the fixed values of O (and the 
other) abundances. Then, in order to check the effect of 
the fixed parameters on the abundance measurement, we 
varied the values of O, Ne, and Fe abundances of the high- 
kT e component from 0.1 to 20 and calculated each best-fit 
parameter. Figure 5 shows the best-fit parameters of Mg, 
Si, S, and Ar abundances for various values as a function 
of O, Ne, and Fe abundance. We also overlaid the x 2 
value on it. We found that the metal abundances depend 
on the O abundance. However, we found the abundances 
for Ar, Si and S are above the solar value. In partic- 
ular, the Ar shows the highest abundance among them. 
Therefore, we concluded that the Ar abundance of the 
Cygnus Loop is significantly higher than the solar value. 
We also confirmed this result from the EPIC spectrum 
(see figure 4 right and table 5). The EPIC result also 
shows a significant high abundance of Ar (SA^l'^) that * s 
consistent with the XIS result within a statistical error. 
Thus, we concluded that the Ar Hea line is detected from 
the Cygnus Loop's ejecta for the first time. 
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4. Discussion 



We investigate the spatial distribution of Ar using the 
Suzaku XIS data. Due to the lack of statistics, we are not 
able to subdivide the XMM-Newton data. We divide the 
XIS FOV into three regions as shown in figure 1. The re- 
gions surrounded by red, blue, and black lines are labeled 
as Region- A, Region-B, and Region-C, respectively. 

Firstly, we applied the bremsstrahlung plus seven 
Gaussian components for each spectrum. The fitting 
method is the same as the one explained in the previ- 
ous section. The results and the best-fit parameters are 
shown in figure 6 and table 6, respectively. From figure 
6, the Gaussian-like structures are seen at about 3.0 keV 
particularly in the Region-A spectrum. We note that the 
line-center energy corresponds to that of the Ar Hea line 
for each fit. The contribution of Ar Hea line emission 
is decreasing from Region-A to Region-C, i.e., from the 
Cygnus Loop center to the rim, suggesting the observed 
Ar is concentrated on the center of the Loop. Figure 7 
shows the radial distributions of the equivalent widths 
for some lines (Mg He/3, Si Hea, S Hea, and Ar Hea). 
From figure 7, the equivalent widths of Si and S lines are 
clearly decreasing from the center to the rim. On the 
other hand, the equivalent width of Mg line is almost con- 
stant in our FOV. According to the previous study (e.g., 
Uchida et al. 2009a), the Si and the S ejecta are concen- 
trated on the Cygnus Loop's center, while the Mg ejecta 
are almost uniform from the center to the rim. Therefore, 
the equivalent- width distributions of Si, S and Mg reflect 
their radial profiles of the ejecta. While the equivalent 
widths of the Ar line have relatively large errors, they are 
also decreasing toward the rim as well as those of the Si 
and S distributions. This result strongly suggests that the 
Ar emission is also derived from the progenitor star. 

In order to confirm it, we next measured the Ar abun- 
dance of each region by applying the two-component 
VPSHOCK model as described in the previous section. 
The spectra from Region-A and Region-B are the super- 
position of both ISM and ejecta components. On the other 
hand, the spectrum from Region-C is virtually dominated 
by the emission from the ISM component. Therefore, we 
applied one-component VPSHOCK model for the Region- 
C spectrum. Figure 8 shows the spectra obtained from 
Region-A (top), Region-B (middle), and Region-C (bot- 
tom). The best-fit parameters are indicated in table 7. 
Each x 2 value shows that two/one-component VPSHOCK 
models are statistically acceptable for Region- A/B and 
Region-C, respectively. From table 7, both the spectra 
from Region-A and Region-B show that the Ar abun- 
dances are significantly higher than that of the solar value. 
Furthermore, the best-fit value of the Ar abundance is de- 
creasing toward the rim as well as those of Si and S. On the 
other hand, the distribution of Mg abundance is almost 
flat. A previous spatially-resolved analysis showed the 
similar distributions of Si and Mg (Uchida et al. 2009a). 
They concluded that these elements are derived from the 
ejecta blown off by a core-collapse SN explosion. From ta- 
ble 7, the co-existence of Si, S and Ar is consistent with a 
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Fig. 6. The 1.3-5.0 keV spectra of the Cygnus Loop. 
Spectral extraction regions are shown in figure 1 top. They 
are fitted with bremsstrahlung plus seven Gaussian compo- 
nents. The residuals are shown in the lower panels. The 
best-fit parameters are shown in table 6. 
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Fig. 7. Equivalent width distributions of Mg He/3 (triangle), 
Si Hea: (X), S Hea (square), and Ar Hea (circle). The hor- 
izontal bars in the x-axis represent the range of each region. 
The upper limit is obtained from the result of Region-C. 

theoretical picture that Si, S, and Ar are mainly produced 
in the static O-burning layer during the stellar evolution 
and the incomplete Si-burning layer during the SN explo- 
sion. Thus we confirmed that the observed Ar emission 
is originated from the ejecta of the Cygnus Loop rather 
than the ISM. Due to the lack of the statistics, we were not 
able to subdivide our FOV to investigate the spatial dis- 
tribution of Ar in more detail. However, considering the 
co-existence of Si and Ar abundances, the Ar-rich ejecta 
may also have a central concentration similar to that of 
Si. 

Some previous studies estimated the progenitor mass 
of the Cygnus Loop to be ~15M (Levenson et al. 
1998; Tsunemi et al. 2007), hence the core-collapse ori- 
gin. However, our results indicate that the O, Ne, and 
Mg abundances are significantly lower than those of Si, 
S, and Ar, which is inconsistent with the standard the- 
oretical core-collapse model. Although there may exist 
O-Ne-Mg-rich region outside our FOV as suggested by 
Uchida et al. (2009a), the low abundances of O-Ne-Mg 
are still remains an open question. Since Ar is produced 
nearer the center of the Fe core than the other lighter ele- 
ments, the abundance measurement will be useful for the 
future estimation of the progenitor mass by expanding the 
observing region and improving the statistics. Figure 9 
shows an example mass fractions of 24 Mg, 28 Si, 36 Ar, and 
56 Fe in a 15M solar-metallicity star as a function of the 
presupernova radius which takes into account explosive 
nucleosynthesis and radioactive decay (Umeda & Nomoto 
2005; Nomoto et al. 2006; Tominaga et al. 2007). The dot- 
ted vertical line represents a boundary between the ejecta 
and the central remnant to yield a canonical amount of 
56 Ni (~O.O8M ; e.g., SN1987A; Blinnikov et al. 2000), 
which locates at R = 1600 km in the progenitor star cor- 
responding to 1.5OM . Since the Ar-rich layer exists near 
the boundary (R > 2080 km), the mass of the ejected Ar 
highly depends on the position of the boundary as well 
as that of Fe. For example, if the mass of the central 
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Fig. 8. The 1.3-5.0 keV spectra of the Cygnus Loop with 
two-component VPSHOCK model (top panels) and one-com- 
ponent VPSHOCK model (bottom panel). The residuals are 
shown in the lower panels. The best-fit parameters are shown 
in table 7. 
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remnant increases only 2% (i.e., 1.53M ), the boundary 
reaches the Ar-rich layer, which may have a large influ- 
ence on the abundance of the ejected Ar. While Si and 
S are also produced near the boundary, these elements 
are abundant as well even in the O-Ne-Mg layer (see fig- 
ure 9). Consequently, the total amount of the ejected Ar 
more significantly depends on the boundary, i.e., progen- 
itor mass than those of Si and S. In general, the Fe abun- 
dance provides us with an important clue to exploring the 
type of the SN explosion. However, as for the Cygnus 
Loop, the low plasma temperature makes it difficult to 
detect the Fe-K line emissions in X-ray. Therefore, the 
accurate measurement of Ar abundance will be important 
in determining the mass of the progenitor star and the 
central compact object that still has not been discovered. 

5. Summary 

In summary, we detected Ar Hea and S He/3 line emis- 
sions from the Cygnus Loop for the first time by accumu- 
lating all the available data of Suzaku and XMM-Newton. 
The Ar abundances obtained from the XIS and EPIC 
spectra are 9.0^3 g and 8.4^2*7? respectively. The high 
metal abundances clearly show that the Ar line originates 
from the eject a of the Cygnus Loop. We also found that 
the central concentration of Ar is closer to the distribution 
of Si and S rather than that of Mg. The measurement of 
Ar abundance will be important in determining the pro- 
genitor mass of the Loop for the future analysis. 
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Fig. 9. Mass fractions of 24 Mg, 28 Si, 36 Ar, and 56 Fe in a 
15 Mq star after taking into account the radioactive decays 
(based on the results of Umeda & Nomoto 2005; Nomoto 
et al. 2006; Tominaga et al. 2007). The dotted vertical line 
represents a boundary between the eject a and the central rem- 
nant of 1.5OM . 
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Table 1. Summary of the 32 Suzaku observations 



Suzaku Observations 



Obs. ID 


Obs. Date 


R.A., Dec. (J2000) 


Effective Exposure 


500020010 


2005-Nov-23 


20 h 56 m 48 s .9, 


31°56 / 54 // .8 


20.4 ks 


500021010 


2005-Nov-24 


20 h 55 m 56 s .O, 


31°56 / 53 // .2 


21.4 ks 


500022010 


2005-Nov-29 


20 h 55 m 05 s .6, 


32°10 / 35 // .4 


21.7 ks 


500023010 


2005-Nov-30 


20 h 54 m 03 s .8, 


32°21 / 47 // .9 


25.3 ks 


501029010 


2006-May-09 


20 h 55 m 00 s .0, 


31°15 / 46 // .8 


13.2 ks 


501030010 


2006-May-10 


20 h 53 m 59 s .3, 


31°03'39".6 


13.9 ks 


501031010 


2006-May-12 


20 h 52 m 58 s .8, 


30°51 / 32 // .4 


18.2 ks 


501028010 


2006-May-13 


20 h 55 m 56 s .3, 


31°28 / 56 // .2 


4.9 ks 


501033010 


2006-May-22 


20 h 50 m 58 s .8, 


30°27 / 00 // .0 


20.0 ks 


501034010 


2006-May-22 


20 h 48 m 49 s .7, 


30°00 / 21 // .6 


13.9 ks 


501032010 


2006-May-25 


20 h 51 m 58 s .6, 


30°39 / 10 // .8 


17.4 ks 


501035010 


2006-Dec-18 


20 h 48 m 16 s .2, 


29°42 / 07 // .2 


11.2 ks 


501036010 


2006-Dec-18 


20 h 47 m 17 s .3, 


30°04 / 21 // .4 


11.8 ks 


501017010 


2007-Nov-ll 


20 h 49 m ll s .3, 


30°59 / 27 // .6 


28.7 ks 


501018010 


2007-Nov-12 


20 h 48 m 18 s .7, 


30°46 / 33 // .6 


21.0 ks 


501019010 


2007-Nov-12 


20 h 47 m 14 s .2, 


30°36'10".8 


16.2 ks 


501020010 


2007-Nov-13 


20 h 46 m 20 s .8, 


30°23 / 22 // .6 


14.7 ks 


501012010 


2007-Nov-14 


20 h 54 m 07 s .6, 


31°57 / 22 // .0 


9.8 ks 


501013010 


2007-Nov-14 


20 h 53 m 08 s .5, 


31°45 / 40 // .3 


16.4 ks 


501014010 


2007-Nov-14 


20 h 52 m 09 s .9, 


31°36 / 43 // .4 


16.9 ks 


501015010 


2007-Nov-14 


20 h 51 m ll s .8, 


31°22 / 08 // .4 


18.3 ks 


501016010 


2007-Nov-15 


20 h 50 m ll s .3, 


31°10 / 48 // .0 


19.3 ks 


503055010 


2008-May-09 


20 h 49 m 48 s .7, 


31°30 / 18 // .0 


22.2 ks 


503056010 


2008-May-10 


20 h 48 m 00 s .0, 


31°10 / 30 // .0 


22.5 ks 


503062010 


2008-May-13 


20 h 56 m 26 s .5, 


30°19 / 55 // .2 


16.9 ks 


503063010 


2008-May-13 


20 h 55 m 16 s .3, 


30°01'44".0 


22.8 ks 


503064010 


2008-May-14 


20 h 53 m 51 s .6, 


29°54 / 42 // .5 


18.2 ks 


503057010 


2008-Jun-02 


20 h 52 m 43 s .8, 


32°26 / 19 // .0 


16.2 ks 


503058010 


2008-Jun-03 


20 h 51 m 17 s .2, 


32°25 / 24 // .6 


19.3 ks 


503059010 


2008-Jun-03 


20 h 49 m 50 s .6, 


32°21 / 50 // .8 


19.5 ks 


503060010 


2008-Jun-04 


20 h 48 m 28 s .2, 


32°17 / 44 // .5 


18.5 ks 


503061010 


2008-Jun-04 


20 h 47 m 22 s .7, 


32°10 / 22 // .8 


26.0 ks 



Table 2. Summary of the 9 XMM observations 



XMM-Newton Observations 



Obs. ID Obs. Date R.A., Dec. (J2000) Effective Exposure 



0082540101 


2002-Nov-25 


20 1 


1 55 m 23 s .6, 


31°46'17" 


.0 


14.7 ks 


0082540201 


2002-Dec-03 


20* 


1 54 m 07 s .2, 


31°30 / 51 // 


.1 


14.4 ks 


0082540301 


2002-Dec-05 


20* 


1 52 m 51 s .l, 


31°15'25" 


.7 


11.6 ks 


0082540401 


2002-Dec-07 


20* 


1 51 m 34 s .7, 


31°00'00" 


.0 


4.9 ks 


0082540501 


2002-Dec-09 


20 l 


1 50 m 18 s .4, 


30°44 / 34 // 


.3 


12.6 ks 


0082540601 


2002-Dec-ll 


20* 


1 49 m 02 s .O, 


30°29'08" 


.6 


11.5 ks 


0082540701 


2002-Dec-13 


20* 


1 47 m 45 s .8, 


30° 13 / 42 // 


.9 


13.7 ks 


0405490101 


2006-May-13 


20* 


1 50 m 32 s .2, 


somw 


.0 


6.5 ks 


0405490201 


2006-May-13 


20* 


1 49 m 54 s .2, 


29°42 / 25 // 


.0 


3.6 ks 
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Table 3. Spectral parameters inferred from the combined XIS and EPIC data (see figure 3) 





XIS (1.3-5.0 keV) 


SlirLKj ^l.O-O.U Ke V J 


Parameter 


Bremss + 6Gauss 


Bremss + 8Gauss 


Rrpmss -1- 7Crrii]ss 


kT ricpvi 

Line Center [eV]: 


61 + - 03 

U ' Oi -0.02 


2.6 (fixed) 

n cq+0.03 




Mg Rea 


1341±1 


1341±1 




Mg He/3 


1572±| 


1572±| 


1555±| 


Si Hea 


1848±1 


1848±1 


1854±1 


Si Re/3 


2178±? 


2176^0 


2185±| 


Si He 7 


2332^ 


2324±20 


2333±21 


S Hea 


2440^4 


2440^4 


2445^ 


S He/3 




2936±g 


2920^° 


Ar Hea 




3128±39 


3132±^ 


X 2 /olof 


1050/998 


993/994 


1214/932 



The errors are in the range Ax 2 < 2.7. 



Table 4. Abundances for \ow-kT e component 



Element Abundance (in units of solar) 



c 


0.46 


N 


0.28 





0.20 


Ne 


0.36 


Mg . . . . 


0.23 


Si 


0.29 


S 


0.18 


Fe 


0.25 


Table 5. Spectral parameters (see figure 4) 


Parameter 


XIS EPIC 


N H [10 2U cm" 2 ] 


2.3 (fixed) 


High-/cT e component: 




kT e [keV] 


°- 69 io.u2 0.49±0.02 


Abundance: 




(=C=N) 


1.42 (fixed) 


Ne 


1.25 (fixed) 


Mg 


0.82±0.03 1.05 (fixed) 


Si 


oi 4 +0.09 4 r 7 +0.12 
6 - i4 -0.10 4 ' ' -0.07 


S 


6 4+ - 5 fi 2 +0 - 4 

°*^-0.4 °' Z -0.2 


Ar 


9.0±|"g 8.4t 2 2 - 5 7 


Fe (=Ni) 


3.20 (fixed) 


Tlower [s CUl" 3 ] 


(fixed) 


log(r upper [s cm" 3 ]) 


10.73±0.04 <12 


Low-kT e component: 




kT e [keV] 


0.28 (fixed) 


Abundances 


fixed* 


ri ower [s cm- 3 ] 


(fixed) 


\og(r upper [s cm" 3 ]) 


11.45 (fixed) 


X 2 /doi 


1096/1004 1020/912 



* Fixed to the values shown in table 4. 
Other elements are fixed to solar values 
The errors are in the range A% 2 < 2.7. 
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Table 6. Spectral parameters inferred from the combined XIS data (see figure 6) 



Parameter Region- A Region-B Region- C 

N H [10 2U cm- 2 ] . . 2.3 (fixed) 

kT e [keV] 0.46±0.01 0.55±0.01 0.35±0.01 

Line Center [eV]: 

MgHea 1339±1 1341^ 1340±? 

MgHe/3 1572±i| 1580^ 1586±? 2 

SiHea 1848±1 1849±1 1848±1 

Si He/3 2172^ 2170t^ 2160±^ 

Si He 7 2304±^ 2332±|| 2348^ 

S Hea 2440±3 2438t^ 2444±^ 

S He/3 2918±|? 2909±i$ 2970±85 

Ar Rea 3104tg 3160±?2 

X 2 /dof 899/957 851/957 1214/932 



The errors are in the range A\ 2 < 2.7. 



Table 7. Spectral parameters inferred from combined XIS data (see figure 8) 





Region-A Region-B 


Region- C 


Parameter 


two-component VPSHOCK one- 


■component VPSHOCK 


N H [10 2U cm-^] 


2.3 (fixed) 




High-/cT e component: 

kT e [keV] 

Abundance: 


0.80^6 0.65±0.01 


U.OO_q q 4 


C 


linked to O abundance 


0.79 (fixed) 


N 


linked to O abundance 


0.20 (fixed) 


O 


1.38 (fixed) 1.46 (fixed) 


0.19 (fixed) 



Ne 

Mg 

Si 

S 

Ar 

Fe (=Ni) 



ric 



s cm 



upper 



[s cm 3 ]) 
Low-kT e component: 

kT e [keV] 

Abundances 

-31 



1.15 (fixed) 

ro+0.24 

i - oz -o.n 
7 2Q+ - 17 

8.96±8;|? 
5.06 (fixed) 



"■3918:85 



1.36 (fixed) 

I. 69±0.02 

9 «Q+0.04 
Z ' O,3 -0.05 
o 94+0.24 

4-lt 2 2 -8 
2.44 (fixed) 

(fixed) 

II. 68±0.03 



0.43 (fixed) 

o-37±5:S 

1 40+ - 81 
<0.33* 
0.42 (fixed) 



11 11+ - 29 

ii - i - L -0.15 



Tlower |S Cm 



log(T M pp er 

X 2 /dof 



s cm 



3 ]) 



0.28 (fixed) 
fixed* 
(fixed) 
11.45 (fixed) 
550/980 565/973 



583/952 



* Fixed to the values shown in table 4. 
''"Estimated by fixing all other parameters. 
Other elements are fixed to solar values. 
The errors are in the range A% 2 < 2.7. 
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